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Abstract 

The goal of the DIRAC experiment at CERN is the lifetime measurement of 
pionium (ir + ir~ atom). Its lifetime is mainly defined by the charge-exchange process 
7r + 7r~ — > 7T°7r . Value of the lifetime in the ground state is predicted in the 
framework of the Chiral Perturbation Theory with high precision: t\s = (2.9 ± 
0.1) • 10~ 15 s. The method used by DIRAC is based on analysis of 7r + 7r _ -pairs 
spectra with small relative momenta in their center of mass system in order to 
find out signal from pionium ionization (break-up) in the target. Pioniums are 
produced in proton-nuclei collisions and have relativistic velocities (7 > 10). For 
fixed values of the pionium momentum and the target thickness the probability 
of pionium ionization in the target depends on its lifetime in a unique way, thus 
the pionium lifetime can be deduced from the experimentally defined probability 
of pionium ionization. Based on ionization cross-sections of pionium with target 
atoms we perform the first direct calculation of the pionium ionization probability 
in the target. 

Introduction 

Pionium is the hydrogen-like Coulomb bound system of two oppositely charged pions. 
Its lifetime is determined by strong 7r + 7r~ — > 7r 7r° annihilation, but it interacts mainly 
electromagnetically with the target atoms, while propagating through the target material 
with the relativistic velocity. The DIRAC experiment at CERN pQ detects 7r + 7r~-pairs 
with low relative momenta in their center of mass system. From the low-momentum part 
of the spectrum, the probability of pionium ionization (break-up) in the target, which 
is the probability for the pionium to be converted into an unbound tt + tt~ pair on the 
exit of the target, is determined. If the dependence of the pionium ionization probability 
in the target as a function of its lifetime is established, this lifetime measurement may 
then be confronted with the value predicted in the framework of the Chiral Perturbation 
Theory [2J. The measurement by DIRAC will test the standard picture of spontaneous 
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chiral symmetry breaking in QCD, providing the experimental evidence in favour of or 
against the existence of a large quark- ant iquark condensate in the QCD vacuum [3]. 

The main thrust of the DIRAC technique is that the dependence of the pionium 
ionization probability in the target as a function of its lifetime can be established with 
sufficient precision. Despite the progress in calculations of ionization and excitation cross- 
sections from different pionium states pE], it proved to be difficult to solve the infinite 
system of differential (kinetic) equations, which describes the evolution of a pionium 
propagating through the target. In the original studies [3 [6] probabilities for the pionium 
to leave the target in a bound state or annihilate were calculated by taking into account 
only low-lying atomic states. The contribution of highly-excited states was estimated 
by extrapolation and the sought-after probability of pionium ionization in the target 
was obtained indirectly by means of unitarity. The precision of the indirect calculation 
was challenged by the fact, that pionium excitation to ever higher lying bound states 
constitutes a major branch in the evolution of pionium. 

In this work we solve the system of kinetic equations by taking into account dynamics 
of highly excited pionium bound states. Within our approach the pionium ionization 
probability in the target is calculated directly, based on ionization cross-sections. The 
formalism of pionium dynamics, described by a set of the probabilistic kinetic equations, 
is reminded in Section [TJ Section [5] is devoted to the direct calculation of the ionization 
probability. 



1 Dynamics of a pionium in the target 

The pionic atoms can be created in inelastic proton-nuclei collisions with the probability 
given by [7j 



(2tt)'|¥(0)| 



dP M dp+dp. 
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where — — fzr is the double inclusive cross section for production of 7r + 7r~ pairs without 
dp + dp_ 

interaction in the final state with both pions produced either directly in hadronic processes 
or through short-lived resonances. P, E and M are the pionium momentum, energy and 
mass in the laboratory frame, respectively, p_ and p + are the momenta of the charged 
pions. Production of the pionium atoms with the angular momentum I > is suppressed. 
The strong interaction between the two pions, which forms the atom, significantly modifies 
|\l/ n o(0)| in comparison to the pure Coulomb wave functions pj. However, only ratios 
between the squares of the wave functions modulus at the origin enter to the following 
calculations. For these ratios the modification due to the strong interaction is almost 
cancelled. The probability of pionium production depends on the principal quantum 
number n as P 
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where \&^ (0) is the pure Coulomb wave function of the 7r + 7r~ atom at zero distance 
and the correction factor (1 + A„) takes into account the effect of finite-size of the pion 
production region and the two-pion strong interaction in the final state. It was found 
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that this correction shifts the probability of ionization on per mille level, therefore we will 
use pure Coulomb wave functions hereafter. For them only nS states are non-zero at the 
origin. If we normalize probability of atom production to unity then the probability of 
atom production in the \nlm) state reads 



^nzJO) 6 l0 5, 



10 



Er=il^ s (o)l 2 ™ 3 C°=iiA 3 ™ 3 C(3)' 

Pioo(O) = 0.832, p 200 (0) = 0.104, p 3 oo(0) = 0.031. 

After production atom can either annihilate (mainly through — > 7r°7r° process^) 
or electromagnetically interact with target atoms. 

The partial decay width of the pionium in IS state into 2tt° is [2]: 



r 2 .o = \o?Jml + - m\ a - Jm»° - ag) 2 m^(l + 8 T ), 6 r = (5.8 ± 1.2) x 10~ 2 . 

The (a° — Oq) difference of the pion-pion S-wave scattering lengths with isospin and 2 
has been calculated [TJ] within the framework of the standard chiral perturbation theory 
(ChPT) [El 

a° -a 2 = (0.265 ± 0.004)m^. 
This difference leads to the predicted value of the pionium lifetime in the ground state 

Tig = (2.9 ±0.1) • 10~ 15 s. (2) 

Lifetime of the atom is inversely proportional to the square of its wave function modulus 

• • mn ■ i l^is(0)| 2 |*fo 

at the origin 141, so m nb states it reads r n s = Ti^rr — — -—- ~ rig r 

Therefore the probability for a pionium with the laboratory momentum P to annihilate 
per unit length is 

f 1 2m^ 

u 1 — 7, = ~ , m nb states, 

nlm I in other states. 

While crossing the target a pionium electromagnetically interacts with target atoms. 
As a result, 7r + 7r _ atom can be either ionized or transit from the initial bound state 
\riilimi) to another bound state |n///m/) (excitation/de-excitation). Hereafter we will 
denote initial and final bound states as \i) and |/), respectively. 

The probability of ionization per unit length from the state \i) is given by 

u/ ion 1 P N A ^ion 

A ion A t ' 

where p is the target density, A is its atomic weight, Na is the Avogadro constant and 
a\ on is the ionization cross section. 

The probability of a pionium excitation per unit length from the state \i) to the final 
state \ f) is given by 

W f = — = 

1 \f A ' 



Another annihilation channel ir + n — > 2-f amounts only about 0.3% [TU1 [TTj . 
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where aj is the discrete (bound-bound) transition cross section. 

The total inelastic cross section gives the probability of an atom to undergo an inelastic 
electromagnetic interaction 

_tot \ _/ i _ion /o\ 
f 

Total cross sections can be calculated owing to the completeness of the eigenstates of the 
Coulomb Hamiltonian. 

Total and transition cross sections for any bound states were initially calculated in 
the Born approximation with the static potential of target atoms [5]. Later a more 
accurate set of cross sections was derived which takes into account relativistic effects, 
multiple photon exchange and target excitations [I] . Moreover, in the latter work authors 
calculated ionization cross sections for any initial bound state with n ^ 8 which provides 
the possibility to perform direct calculation of the pionium ionization probability in the 
target, which is the subject of the present work. From general formulas for transition 
cross sections in both sets [HI Hj one can derive that for any states \i) and |/) o{ = a\. 
Comparison between different sets of cross sections was performed [6], where authors 
found that uncertainties in most precise sets of cross-sections for Ni target will cause only 
1% uncertainty in the pionium lifetime. Uncertainty due to the accuracy of cross sections 
is expected to dominate precision of the ionization probability dependence on the pionium 
lifetime. 

The dynamics of the pionium interaction with target atoms is supposed to be described 
by a set of kinetic equations [5] using the probabilities p%(s) to find the n + ir~ atom in 
the definite quantum state \i) at a distance s from the production point. This approach 
ignores any interference between different pionium bound states. For low n most of 
interference effects are suppressed at typical pionium momenta in DIRAC (3 8 GeV/c) 
as the mean free path between pionium inelastic interactions is usually longer than the 
formation time of atomic system multiplied by its velocity. Nevertheless even for low n 
some interference effects can take place due to the accidental degeneracy of energy levels 
of hydrogen-like atoms. This problem was considered in the framework of the density 
matrix formalism [15]. It was found that the interference between quantum states with 
small n does not change the result based on a set of probabilistic kinetic equations (their 
difference is less than per mille). 

Eigenstates of the Coulomb Hamiltonian form a countable set of discrete levels. For 
numerical calculations we will take into account only levels with a principal quantum 
number n ^ n max . For a given principal quantum number n there are n 2 states \nlra) 
with different orbital and magnetic quantum numbers. We will denote the total number 
of discrete bound states taken for the calculations as N. To make the system complete 
we introduce a cross section cr" which stands for the sum of transitions from state \i) to 
any discrete state with a principal quantum number rif > n max : 

f:n f >n 

max 

It is straightforward to write the probability of pionium production in all bound states 
above n max : 

ax 
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Finally we will write the system of kinetic equations in the matrix form 
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Diagonal terms describe the total decrease of the level population 

Wf = 



P N A 

^tot Ti/anh 



(5) 
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System (jSJ) is a system of linear homogeneous differential equations of order 1 with 
constant coefficients. The rank of the matrix is N, with three low lines being a linear 
combination of the first iV lines. It is exactly solvable 
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where Ai, . . . , are eigenvalues and aS k > their corresponding eigenvectors. Symmetry of 
the upper left N-by-N corner guaranties that all its eigenvalues are real [TS] . Coefficients 
Cfc are fixed from initial conditions ([I]), (Hj). The probability of pionium ionization at the 
distance s after the production point is expressed through the solution ([6]) 



Pi. 



-1) 



wr n a 



(k) 



DIRAC uses very thin targets (their nuclear efficiency is less than 10 -3 ), therefore 
atoms are produced nearly uniformly over the target thickness sq. Hence the probability 
for a pionium to leave the target in the state \i) reads 



1 f s ° 

Pi{s ) = — Pi(s)ds = 

s o Jo V 



c k a 



(k) 
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while the probability of ionization on the exit of the target is 



[so) 



■ on W 



(7) 



Expressions for the probability of annihilation P ana and for the probability P u to reach 
any excited state with n > n max on the exit of the target have the same form as ((7|) if one 
substitutes W\ on with W? nh and respectively. 

In Tabled] we illustrate this solution as a function of n max for a pionium atom produced 
in 95 jum thick Ni target with the momentum P = 4.6 GeV/c, corresponding to the 
average laboratory momentum of produced pioniums in the kinematic range of the DIRAC 
experiment. Eigenvalues were numerically found by functions from the LAPACK [T7] 
package. System of equations (jSJ) is constructed in a way that as soon as an atom reaches 
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Table 1: Results for different probabilities on the exit of the target as a function of n 
if the dynamics of highly excited states (with n > n max ) is not taken into account 



^max 


pA 

dsc 


pA 


P A n 




5 


0.0928 


0.4407 


0.2438 


0.2227 


6 


0.0937 


0.4407 


0.2586 


0.2070 


7 


0.0943 


0.4407 


0.2715 


0.1935 


8 


0.0947 


0.4407 


0.2828 


0.1817 



a state with n > n max it effectively quits from calculations and is kept intact, though 
in reality it is expected to undergo further electromagnetic interactions, e. g. it can be 
ionized or de-excited to the low-lying states. Therefore, P,£, c is at least the probability 
for a pionium to leave the target in any bound state with the principal quantum number 
n ^ ^max- Pu is the probability for atoms to reach states with n > n max , which amounts 
to about 20%. This numerical value is in agreement with the earlier calculations [6J 
(Fig. 3(d)). Numerical precision of the above solution can be estimated from the inequality 

1 1 — P A — P A — P. A — P A I <r 1 . 1 n -12 

I x ± dsc ± anh ± ion ' u | x 1 xw > 

thus round-off errors do not affect the result. 

As an atom transits to a state nf > Hi its effective radius of electromagnetic interac- 
tions grows and its characteristic ionization length is getting shorter 

\ ion \ ion ~ 9 nm 

A |n>n max ,Zm) ^ A |n ma x=8,Zm) ~ Z 

The target used in DIRAC is 95/im thick, therefore highly excited atoms have a chance to 
leave the target in a bound state only if they were created close to the exit of the target. 
Otherwise these highly excited atoms will be ionized. This allows us to set the range for 
the ionization probability in the target: 

0.2828 = P£ n < P ion < P£ n + P* = 0.4645. (8) 

Here the upper bound corresponds to the case when all highly excited atoms are ionized, 
while the lower bound P { A n is at least probability of ionization from states with n ^ n max . 

From Table [T] one can conclude that above upper and lower bounds converge slowly 
with increase of n max and in this way it would be difficult to increase n max in order to 
achieve precision required by DIRAC (per cent level). 



2 Evolution of highly-excited states 

Rather than trying to solve the system (j^j) directly we will modify it in order to get 
the lower bound of the ionization probability by taking into account dynamics of highly 
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excited states with n > n r 
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(9) 



Here 



pN A 



(J, 



-l.min 



A - n. m ax+l,min) u n max +l,min j/J^/ V l n max+l,i'm') . 

is the lower bound of the probability of ionization per unit length from any state with 
n > n max, because the ionization cross section tends to grow with increasing of the 
principal quantum number n due to the corresponding expansion of the atomic radius. The 
minimal and maximal values of the ionization cross-section for different principal quantum 
numbers are drawn in Fig. [TJ To find the lower bound of the ionization probability, further 
we require all probabilities of ionization per unit length from any state \nlm) do not exceed 

TXrion. 
' ' u ' 

pN A 

; n ma x+l,min. 



mm), mm 



-<J, 



a J '\nlm) ,min' ^ \nlm) , mm L1LL1L \ u \nlm) > 

The diagonal term Wj, which describes the level de-population, is changed accordingly 
to fulfill ©• 

Upper bound of the probability of de-excitation per unit length from all states with 
n > n max to a state |/) with n/ ^ ra max is obtained from the following inequality 

w/ Pi < w !- E Pi = wip u = w? Pu . 

V.Tl\ ^^max I'.Tl^ >?lrriax i:7l^Tl m ax 

Finally the diagonal term, which describes de-population of discrete states with n > n max , 
is 



mm < err. 



where W* nh 



' * u 



2m„ 



pN A 
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anh 



is the upper bound of the probability for a pionium to 



^T-ls(rimax + l) 3 

annihilate from any state with n > n max per unit length. The rank of the new system 
is N+l. The system ([9]) is constructed in the way that ionization is underestimated and 
all competitive processes including de-excitation from high n states (thus transitions to 
bound states with even lower ionization) are overestimated, therefore the solution is the 
mathematical lower bound of the probability of ionization. Numerical results for different 
n mH are presented in Tab. [21 

Upper (jHJ) and lower bounds effectively squeeze the solution (Fig. EJ), for n mSuX = 8 they 

are 

nmax r?min 

,, > (1Q) 



4567 = P. B < p < P. A + P 1 



0.4645, 



IP- 



0.8-10" 



Precision of the above calculated value of the pionium ionization probability in the 
target is comparable to the ~ 1% uncertainties in the P ion value due to precision of 
the corresponding electromagnetic cross-sections [B]. 
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Table 2: Numerical solution of the system ([9]) for the lower bound of the probability of 
pionium ionization in the target function of n max 



^max 


pB 

dsc 


P mih 






5 


0.1109 


0.4416 


0.4468 


0.00067 


6 


0.1068 


0.4411 


0.4517 


0.00030 


7 


0.1041 


0.4409 


0.4548 


0.00015 


8 


0.1023 


0.4408 


0.4567 


0.00008 



Upper and lower bounds of the probability of pionium ionization in the target as 
a function of its lifetime in the ground state are shown in Fig. [3j The corresponding 
uncertainties are shown as dashed lines around the value of the pionium lifetime (j2J) 
predicted by theory. The DIRAC collaboration reported the measured value of P mn = 
0.4521q;q39 [IB] , based on part of the collected data. While further analysis will reduce 
uncertainties of the experimental result, uncertainties of the solution (ITU]) are expected to 
be within precision, required by DIRAC. Range can be further shrunk by extrapolation 
as shown in Fig. El 

We have to emphasize that upper and lower bounds of the ionization probability 
squeeze the solution with required precision due to the fact, that for atoms with a principal 
quantum number n > 8 a characteristic ionization length is less than 2 //m, which is much 
shorter than the target thickness of 95 /im. If one selects very thin target (e.g. 10 /im 
thick Ni) then the upper and lower bounds will show worse convergence (Fig. H]). 

Conclusions 

We derived a mathematical approach to solve a system of kinetic equations, which de- 
scribes evolution of relativistic 7r + 7r~ atoms propagating through the target. In this ap- 
proach we reduce the system of kinetic equations, which formally contains infinite number 
of equations, to the finite set of equations, which is solved exactly. The solution represents 
lower and upper bounds for the probability of pionium ionization in the target. These 
lower and upper bounds effectively squeeze the solution to the value of the probability 
of ionization with 1% precision, which is within requirements of the DIRAC experiment. 
Thus the first direct (based on ionization cross sections) calculation of the probability of 
ionization has been performed. We confirm that the contribution of highly-excited states 
(with the principal quantum number n > 8) to the probability of ionization is significant 
(>l/3). 

Author would like to thank L. Afanasyev, L. Nemenov, A. Tarasov and V. Yazkov for 
many discussions about the problem. 
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Figure 1: The minimal ^rnin {(7|°° m ^}J and the maximal ^max {c[™' m ')} J pionium 

ionization cross sections on the Ni atom for different initial principal quantum numbers n 
of the pionium (according to [1]). 




Figure 2: Upper and lower bounds of the probability of pionium ionization in the target 
Pi on as a function of n max , fitted by ae anmax + c functions to guide the eye (a, a and c are 
free parameters). 
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Figure 3: Upper and lower bounds of the pionium ionization probability in the target Pj 
as a function of the pionium lifetime in the ground state. 
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Figure 4: Upper and lower bounds of the pionium ionization probability in the target P\ 
as a function of n max for 10 /zm thick Ni target. 
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